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Abstract. Spin- 1/2 Ising-Heisenberg model with XYZ Heisenberg pair interaction and two different Ising quartic interactions 
is exactly solved with the help of the generalized star-square transformation, which establishes a precise mapping equivalence 
with the corresponding eight-vertex model on a square lattice generally satisfying Baxter's zero-field (symmetric) condition. 
The investigated model exhibits a remarkable weak-universal critical behavior with two marked wings of critical lines along 
which critical exponents vary continuously with the interaction parameters. Both wings of critical lines merge together at 
a very special quantum critical point of the infinite order, which can be characterized through diverging critical exponents. 
The possibility of observing reentrant phase transitions in a close vicinity of the quantum critical point is related to a relative 
strength of the exchange anisotropy in the XYZ Heisenberg pair interaction. 

Keywords: Ising-Heisenberg model, eight-vertex model, quantum critical point, weak universality, reentrant phase transitions 
PACS: 05.50.+q, 05.70.Fh, 05.70.Jk, 64.60.De, 64.60.F-, 75.30.Gw, 75.30.Kz 

INTRODUCTION 

Exactly soluble quantum spin models belong to the most fascinating topics to deal with in the area of modern 
equilibrium statistical mechanics [110,[3t]. It should be pointed out, however, that quantum effects usually compete 
with a cooperative nature of spontaneous long-range ordering and thus, it is quite intricate to find an exactly solvable 
model that simultaneously exhibits both spontaneous long-range order as well as obvious macroscopic features of 
quantum origin. On the other hand, it is a competition between quantum and cooperative phenomena that is an essential 
ingredient for observing a quite remarkable and unexpected behavior of low-dimensional quantum spin models. 

The hybrid Ising-Heisenberg models on decorated planar lattices, whose nodal sites are occupied by the classical 
Ising spins and decorating sites by the quantum Heisenberg ones, belong to the simplest rigorously solved quantum 
spin models that exhibit a spontaneous long-range ordering with apparent quantum manifestations. It is worthwhile to 
remark, moreover, that the Ising-Heisenberg planar models |,4,. 5, 6,. 7. 8. 9. 10, 11. 12. . 131 where a finite cluster of the 
Heisenberg spins interacts with either two or three nodal Ising spins are in principle tractable by the use of generalized 
decoration-iteration or star-triangle transformations lITil ITsl Il6ll . which establish a precise mapping equivalence 
between them and the spin-1/2 Ising model on the corresponding undecorated planar lattice |1, 15, 17, 18, 19, 20]. 
Among other matters, the exact solutions for this special class of the Ising-Heisenberg planar models might serve in 
evidence that these rigorously solvable models exhibit a strong-universal critical behavior, which can be characterized 
by critical exponents from the standard Ising universality class. Contrary to this, the more interesting weak-universal 
behavior li2lil of the critical exponents has been recently announced for two Ising-Heisenberg planar models i22il23ll . 
where a finite cluster of the Heisenberg spins interacts with four nodal Ising spins. The spin-1/2 Ising-Heisenberg 
model with the pair XYZ Heisenberg interaction and two quartic Ising interactions ll23ll has surprisingly turned out 
to be the fully exactly solvable model due to a validity of the precise mapping equivalence with Baxter's zero-field 
(symmetric) eight-vertex model I ll l24l l25ll . The main purpose of this work is to examine in detail how the weak- 
universal critical behavior of this exactly soluble model depends on a spatial anisotropy in two quartic Ising interactions 
and on the exchange anisotropy of the XYZ Heisenberg pair interaction, whose effect have not been dealt with in our 
preceding work Ii23il for the most general case. 



FIGURE 1. The elementary unit cell of the spin- 1/2 Ising-Heisenberg model. Full (empty) circles denote positions of the 
Heisenberg (Ising) spins, thick solid line represents the pairwise XYZ Heisenberg interaction between the apical Heisenberg spins 
and both types of broken lines connect spins involved in the quartic Ising interactions. Thin solid lines connecting four Ising spins 
are guide for eyes only. 

This paper is so organized. In the following section, we will describe the hybrid Ising-Heisenberg model and recall 
basic steps of the exact mapping procedure to the zero-field eight-vertex model. The most interesting results for the 
ground-state and finite-temperature phase diagrams, which are supported by a detailed analysis of critical exponents, 
are subsequently presented in the next section. Finally, some concluding remarks are mentioned along with a brief 
summary of the most important scientific achievements in the last section. 

ISING-HEISENBERG MODEL AND ITS EXACT MAPPING EQUIVALENCE WITH 
THE ZERO-FIELD EIGHT- VERTEX MODEL 

Let us consider a two-dimensional lattice of edge-sharing octahedrons. Figure[T]schematicaIIy displays the elementary 
unit cell of the two-dimensional lattice, i.e. an octahedron, which contains four Ising spins a = 1/2 in its basal plane 
and two Heisenberg spins 5 = 1/2 in its apical positions. Let each edge of the octahedron, which connects two Ising 
spins, be a common edge of two adjacent octahedrons so that an ensemble of all Ising spins will form a square 
lattice and the Heisenberg spins will be located above and below a center of each elementary square face. Suppose 
furthermore that both apical Heisenberg spins interact together via the pair XYZ Heisenberg interaction, while they 
also take part in two different quartic Ising-type interactions with two Ising spins from opposite corners of a square 
face (see figure [U. For further convenience, the total Hamiltonian can be written as a sum over all elementary unit 
cells (octahedrons) — '-^^p, where each cluster Hamiltonian Jfp contains one pair interaction between the apical 
Heisenberg spins and two quartic Ising-type interactions between the Heisenberg spins and their four Ising neighbors 

Above, the interaction parameters Jx,Jy,J^ label spatial components of the anisotropic XYZ interaction between the 
Heisenberg spins, while the interaction parameters Ji and J2 label two quartic Ising-type interactions between both 
apical Heisenberg spins and two Ising spins from opposite corners of a square face along two different diagonal 
directions (see figure[T]i. 

The crucial step of our calculation lies in an evaluation of the partition function. A validity of the commutation 
relation [Jfj, — between different octahedron-cluster Hamiltonians allows a partial factorization of the partition 
function into the following product 

^iHM = E nTi>exp(-j34,) = E l\^picT;„a;2,CT;^,CT;^, (2) 

{a} P {a} P 

where j3 = 1/ (kgT), is Boltzmann's constant and T is the absolute temperature. The summation L{c7} to emerge 
in Eq. (|2|i is carried out over all possible configurations of the Ising spins, the product runs over all octahedron unit 
cells and the symbol Trp stands for a trace over spin degrees of freedom of the Heisenberg spin pair from the pth 
octahedron. In the latter step of our calculation we have used a straightforward diagonalization of the Hamiltonian ([TJ 
of the pth octahedron in order to obtain the relevant trace over spin degrees of freedom of the Heisenberg spin pair. 



This procedure yields some effective Boltzmann's weight cOp, which explicitly depends merely on four Ising spins Opi, 
<7p2, CpS and ffp4 from the basal plane of the pth octahedron. In addition, the explicit form of the effective Boltzmann's 
factor (Op immediately implies a possibility of performing the generalized star-square transformation 
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which substitutes the effective Boltzmann's weight (Op by the equivalent expression containing two pair (Ri and R2) 
and one quartic {R4) interaction between four nodal Ising spins from an elementary square face of the square lattice. 
Of course, the algebraic transformation Q must satisfy the 'self-consistency' condition, which means that it must hold 
independently of spin states of four Ising spins involved therein. It can be easily verified that a substitution of all sixteen 
possible spin configurations of the nodal Ising spins gives just four independent equations, which unambiguously 
determine so far not specified mapping parameters Rq, Ri, R2, and R4 
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that can be expressed in terms of four different Boltzmann's weights COi{i = 1,3,5,7) 
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Note that the effective Boltzmann's weights ([Sj-® are assigned to eight spin configurations of four nodal Ising spins 
explicitly specified in round brackets (the symbol ± denotes the spin state ct^q, = ±1/2, a = 1,2,3,4), as well as, 
another eight spin configurations, which can be obtained from them under the reversal of all four nodal Ising spins. 

A substitution of the generalized star-square transformation (|3]l with appropriately chosen mapping parameters (|4|i 
into Eq. (|2]i then straightforwardly leads to a precise mapping equivalence between the spin- 1/2 Ising-Heisenberg 
model and the zero-field eight-vertex model. As a result of this procedure, one actually obtains a simple mapping 
relationship 
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which connects the partition function of the investigated Ising-Heisenberg model with the partition function of 
the spin- 1/2 Ising model on two interpenetrating square lattices that are coupled together by means of the quartic 
interaction (A^ is the total number of the Ising spins). It should be mentioned that the latter model is nothing but one 
of many alternative definitions of the zero-field eight-vertex model U H^JEsj], which is reformulated in the Ising 
spin representation following the ideas of Wu |26], Kadanoff and Wegner [27]. Besides, the physical meaning of the 
mapping parameters Ri, R2 and R4 becomes quite evident from the mapping transformation ([3]). The parameters Ri 
and R2 denote the effective pair interactions in two different interpenetrating Ising square lattices and the mapping 
parameter R4 determines the effective quartic interaction that couples together both Ising square lattices. Last but 
not least, the mapping parameter Rq (multiplicative factor in Eq. ^) in fact represents the partition function of the 
Heisenberg spin pair in the effective field produced by the four enclosing nodal Ising spins. 

It becomes quite clear from the mapping relationship (|9]l that the Ising-Heisenberg model becomes critical just if the 
corresponding zero-field eight-vertex model becomes critical as well. As a matter of fact, the partition function SIum 
will exhibit a non-analyticity if and only if the corresponding partition function ifg-i' will exhibit the same kind of 



the non-analytic behavior, because the parameter Rq is analytic in the whole region of interaction parameters. In this 
regard, the respective critical points of the Ising-Heisenberg model can be found from the relevant critical condition of 
the zero-field eight- vertex model 111 124112511 

0)1 + a>3 + 0)5 + to? =2max{ 0)1,0)3, 0)5, tOy}, (10) 

which determines phase transitions of the Ising-Heisenberg model on assumption that the effective Boltzmann's 
weights ©-dD are substituted into this critical condition. It should be also mentioned that the critical exponents 
of the zero-field eight- vertex model IH llftl continuously change with the parameter jX—2 arctan( 0)5 0)7/ 0)i 0)3 ) ' 
by following the formulas 

H I61J. 2jx 8/x 4 

Consequently, the relations (fTTl i will also govern changes of the critical exponents of the Ising-Heisenberg model when 
the effective Boltzmann's weights ©-([Sli are used for a calculation of the parameter ji. 



RESULTS AND DISCUSSION 

In this part, let us proceed to a discussion of the most interesting results obtained for critical properties of the spin- 1/2 
Ising-Heisenberg model with the pair and quartic interactions. Before doing so, however, it is worthy of notice that 
some special cases of this model system have already been investigated by the present authors in our earlier paper 
|j25], where the critical behavior of the particular case with two identical quartic interactions J\ — J2 and the more 
symmetric XXZ Heisenberg interaction {Jx = Jy^ J^) was explored in detail by assuming both ferromagnetic as well 
as antiferromagnetic pair interaction. Exact results for this more symmetric version of the Ising-Heisenberg model 
indicate that the model with the antiferromagnetic pair interaction exhibits less significant changes of both critical 
temperatures as well as critical exponents than the model with the ferromagnetic pair interaction. As a matter of fact, it 
has been demonstrated that only the Ising-Heisenberg model with the ferromagnetic pair interaction shows a unusual 
quantum critical point of the infinite order, which characterizes a remarkable singular behavior of the critical exponents 
near the isotropic limit of the Heisenberg pair interaction (J^ = Jy = J^). With this background, the main purpose of the 
present work is to shed light on how this strange weak-universal critical behavior will change by introducing a spatial 
anisotropy in two quartic Ising-type interactions (Ji 7^ J2) or by assuming the most general form of the XYZ exchange 
anisotropy in the Heisenberg pair interaction. 

It is worthwhile to remark that the investigated Ising-Heisenberg model possesses a rather high symmetry, because 
all four different Boltzmann's weights (|5]l-([8]) are mutually interchangeable under the transformations Ji —Ji and/or 
J2 ~J2- This means, among other matters, that one may further restrict both quartic interaction parameters to positive 
values, since the transformations J\ —Ji and J2 —J2 merely cause rather trivial changes of the nodal Ising spins 
((^Pi.(^p2'(^p3'(^p4) ^ (±(^pi.(^p2.TCT;3'f^p4) and {ol^, 01^,01^,0;^) ^ (ct^i , ia^j, (T^g , T(t;4), respectively. Owing 
to this fact, let us further assume that all the interaction terms J^, Jy, Jz, Ji, and J2 entering the Hamiltonian ([T]i are 
positive and moreover, the symmetry of the Hamiltonian allows us to consider > Jy and Ji > J2 without loss 
of the generality. For easy reference, the number of free parameters is lowered by introducing the following set of 
dimensionless parameters: k^T /J, marks the dimensionless temperature, Ji /J, and J2 /Jz are proportional to a relative 
strength of the quartic Ising-type interactions, and finally, the parameters Jx/Jz and Jy/Jx measure a relative strength 
of the exchange anisotropy in the XYZ Heisenberg pair interaction. The former anisotropy parameter Jx/Jz determines 
a difference in the exchange interactions along a quantization z-axis and its perpendicular x-axis, whereas the latter 
anisotropy parameter determines a difference between the stronger and weaker exchange interaction in the xy-plane. 

First, let us take a closer look at the ground-state behavior. It can be readily understood that the ground-state spin 
arrangement will always correspond to the lowest-energy eigenstate that enters into the greatest Boltzmann's weight 
among the four Boltzmann's weights given by Eqs. (|5ll-([8]l. In the zero temperature limit, the greatest Boltzmann's 
weight is 0)i if Jy < Jz or 0)3 if Jy > Jz- This observation would suggest that the ground-state spin arrangement will 
basically change whenever the weaker exchange interaction Jy in the xy-plane exceeds the exchange interaction 7, 
along the quantization axis. One actually finds that the lowest-energy eigenstate is either 
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if the reverse inequality holds. In Eqs. (fT2l) -(fT3ll. the former ket vector unambiguously determines the states of four 
nodal Ising spins from an elementary square face of the pth octahedron and the latter ket vector unambiguously 
specifies the relevant state of the Heisenberg spin pair. It should be also noticed that another equivalent representations 
of these eigenstates can be obtained from the eigenvectors (fT2]i-(fT3]l under the reversal of all four nodal Ising spins 
and consequently, the phases |I) and |II) are both four-fold degenerate. 

Let us now make a few comments on spin arrangements appearing in both ground-state phases. In the phase |I), 
the Ising spins placed on a square lattice either exhibit a perfect ferromagnetic or antiferromagnetic long-range order, 
which is accompanied with the entangled spin state (|+,+) + | — ))/ \/2 consisting of both ferromagnetic states of 
the Heisenberg spin pairs. On the other hand, the Heisenberg spin pairs reside in the phase |II) the entangled spin state 
(1 + ,— ) + | — ,+))/ \/2 consisting of both antiferromagnetic states and the Ising spins prefer a superantiferromagnetic 
long-range order with the ferromagnetic aUgnment in a horizontal direction and the antiferromagnetic alignment in 
a vertical direction, or vice versa. To compare with, is might quite useful to mention that the phase |II) preserves its 
character even if the less general case of the XXZ exchange anisotropy is assumed, whereas the spin arrangement of 
the phase |I) drastically changes disentanglement of both ferromagnetic states of the Heisenberg spin pairs 
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In the particular limit = Jy, the phase |I') actually becomes macroscopically degenerate as the Heisenberg spin pairs 
may choose independently of each other one of two ferromagnetic states |±, i)^^. Hence, it follows that the phase |I') 
will exhibit a rather high macroscopical degeneracy of the order 4 x 2^ , which is proportional to the total number of 
the Heisenberg spin pairs (remember that simultaneously marks the total number of the Ising spins, the total number 
of the Heisenberg spin pairs, as well as, the total number of elementary unit cells). 

Now, we will turn our attention to a detailed study of finite-temperature phase diagrams. One should recall that phase 
transition lines of the spin- 1/2 Ising-Heisenberg model with the pair and quartic interactions can be straightforwardly 
calculated from the critical condition of the corresponding zero-field eight-vertex model by substituting the effective 
Boltzmann's weights (|5j-([8]l into Eq. ( fTOl i. It is worthwhile to remember, moreover, that the greatest Boltzmann's 
weight is either (£i\ or 03. For both particular cases, the critical condition can uniquely be expressed as 
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where jSc = l/{kBTc), is the critical temperature and the plus (minus) sign applies for a situation when coi > Oh, 
(0)1 < 0)3). 

By exploiting the critical condition ( fTST i, let us examine first how a spatial anisotropy in two quartic Ising-type 
interactions influences the critical behavior of the model under investigation. To see this effect, the critical temperature 
is plotted in figure|2]against the anisotropy parameter for the particular case of the XXZ exchange anisotropy (i.e. 
J.x = Jy) when a strength of the stronger quartic interaction is fixed {Ji/Jz = 1 -0) and a strength of the weaker quartic 
interaction varies. As one can see, the phase diagram consists of two marked wings of critical lines that separate 
both spontaneously long-range ordered phases |I) and |II), since both critical lines merge together at the ground-state 
boundary J^ = Jy = J^ between these two phases. It is noteworthy that the left (right) wing represents a line of critical 
points of the phase |I) (the phase |II)), which can be obtained as a numerical solution of the critical condition ( fTSl l 
by assuming the plus (minus) sign therein. It becomes quite evident from figure |2] that one may also found reentrant 
phase transitions slightly above the ground-state boundary between the phases |I) and |II), i.e. under the constraint 
Jx = Jy > Jz- In this parameter region, the reentrant phase transitions from the spontaneously ordered phase |II) to the 
disordered phase and from the disordered phase to the spontaneously ordered phase |I) take place because the free 
energy of the phase |I) decreases much more rapidly with the increase of temperature than the free energy of the phase 
I II) due to the much higher entropy gain of the phase |I). Another interesting fact to observe here is that a spatial 
anisotropy in the quartic Ising-type interactions does not qualitatively affect the critical behavior of the investigated 
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FIGURE 2. The dimensionless critical temperature as a function of the anisotropy parameter Jx/Jz for the Ising-Heisenberg 
model with the XXZ exchange anisotropy (Jx = A ) when a strength of the one quartic Ising-type interaction is fixed (J\/Jz = 1-0) 
and that of another one quartic Ising-type interaction J2/J: varies. 

model system. As it can be easily understood from Eqs. (|4]i and ©-([HI), the spatial anisotropy in the quartic Ising- 
type interactions Ji ^ J2 merely causes an anisotropy in the effective pair interactions Ri ^ R2 of two Ising square 
lattice coupled together through the effective quartic interaction R4 when speaking in the language of the Ising spin 
representation of the corresponding zero-field eight-vertex model. Of course, it is well known that the anisotropy in 
those pair interactions does not fundamentally affect the overall critical behavior of the zero-field eight-vertex model. 

For the sake of comparison, we have displayed in figure [3] the critical temperature as a function of the exchange 
anisotropy Jx/Jz for one illustrative example of the less symmetric XYZ exchange anisotropy {Jy/Jx — 0.5) when 
a strength of the stronger quartic interaction is fixed {J\/Jz ~ 10) and a strength of the weaker quartic interaction 
varies. It can be clearly seen from this figure that the finite-temperature phase diagram of this more general case quite 
closely resembles the one formerly discussed by the analysis of the particular case with the XXZ exchange anisotropy. 
Actually, one still finds two wings of the critical Unes that separate both spontaneously long-range ordered phases 
|1) and 1 11), whereas the zero-temperature transition between these phases moves towards the higher values of the 
exchange anisotropy Jx/Jz it occurs just if the weaker exchange interaction in the xy-plane overwhelms the one 7, 
along the quantization axis. However, the most obvious difference between the phase diagrams shown in figures|2]and 
[3] is that the latter phase diagram does not imply an existence of the reentrant phase transitions near the ground-state 
boundary between the phases |1) and |11). This specific feature can be attributed to the fact that the macroscopical 
degeneracy of the phase |1) is entirely lifted whenever the most general case of the XYZ exchange anisotropy is 
assumed. Hence, it follows that both wings of critical lines tend to zero temperature with an infinite gradient and thus, 
there does not emerge reentrant phase transitions near the ground-state boundary between the phases |1) and |11). To 
provide an independent check of the aforementioned scenario, the critical temperature is plotted in figure|4]against the 
exchange anisotropy Jx/Jz for one selected value of the quartic Ising-type interactions {J\/Jz = Ji/Jz — 1-0) four 
different strengths of the exchange anisotropy Jy/Jx- Figure |4] apparently supports all previous statements, since the 
critical lines obviously have vertical tangents at the zero-temperature transition between the phases |1) and |11), which 
is generally shifted towards the higher values of the anisotropy parameter Jx/Jz upon strengthening the exchange 
anisotropy in the xy-plane (i.e. when the ratio Jy/Jx is lowered). Altogether, it could be concluded that the spatial 
anisotropy in the quartic Ising-type interactions has a less significant impact on the overall critical behavior compared 
to the exchange anisotropy of the Heisenberg pair interaction. Therefore, our subsequent analysis will be mainly 
concentrated on how the exchange anisotropy in the Heisenberg pair interaction influences the most essential features 
of the critical behavior 

At this place, let us make a few remarks on possible changes of the critical behavior and critical exponents, which 
can be induced upon varying the anisotropy parameters Jx/Jz Jy/Jx- For this purpose, the figures |5}^7] display 
typical changes of the critical temperature and the critical exponent a with the anisotropy parameter Jx/Jz for the 
Ising-Heisenberg model with a fixed strength of the quartic Ising-type interactions {J\/Jz — Ji/Jz— 10) three 
different values of the exchange anisotropy Jy/Jx = 10, 0.9 and 0.5, respectively. In these figures, solid lines scaled 





FIGURE 3. The dimensionless critical temperature as a function of the anisotropy parameter JxlJz for the Ising-Heisenberg 
model with Ihe XYZ exchange anisotropy (J^ 7^ Jy) when a strength of the one quartic Ising-type interaction is fixed {J\/Jz = 1-0) 
and that of another one quartic Ising-type interaction J2IJZ. varies. 




FIGURE 4. The dimensionless critical temperature as a function of the anisotropy parameter Jx/Jz for the spin- 1/2 Ising- 
Heisenberg model with a fixed strength of the quartic Ising-type interactions {J\/Jz = Jl/Jz — 1-0) four different values of 
the exchange anisotropy Jy/Jx- 



with respect to left axes depict a variation of the critical temperature with the exchange anisotropy Jx/Jz, while broken 
lines scaled with respect to right axes show in a semilogarithmic scale the relevant changes of the critical exponent 
a. The figures on the left show the overall finite-temperature phase diagrams, whereas the figures on the right show 
in an enlargened scale the most striking part of these phase diagrams in a close vicinity of the ground-state boundary 
between the phases |1) and |11). 

It is quite obvious from figures|5}^7]that the critical exponent a varies continuously along the line of critical points, 
because its value changes with the interaction parameters involved in the Hamiltonian ([U according to the relation (fTTI) 
through the respective changes of the parameter /i. It should be pointed out, moreover, that the relevant changes of the 
critical exponent bring also insight into a nature of phase transitions as the order of phase transitions is proportional 
to r = 2 — a (see for instance pp. 16-17 in Reference 1 1]). From this point of view, the displayed critical lines turn 
out to be lines of rather smooth continuous phase transitions, since the respective variations of the critical exponent a 
are restricted to the range a E (—00,0) and consequently, the order of phase transitions is r > 2. However, the most 
striking finding to emerge from figures|5]-[7]is that the critical exponent a exhibits a very special singular behavior in 




FIGURE 5. The changes of the critical exponent a along the critical line of the spin- 1/2 Ising-Heisenberg model with a fixed 
relative strength of the quartic Ising-type interactions Ji jj- = J^lh. = 1-0 and the exchange anisotropy JyjJx = 1.0. The solid line, 
which is scaled with respect to the left axis, shows the critical temperature as a function of the exchange anisotropy Jx/J?- The 
broken line, which is scaled with respect to the right axis, displays in a semilogarithmic scale the relevant changes of the critical 
exponent a along this critical line. Figure 5b) shows a detail of the phase diagram near the zero-temperature transition. 




FIGURE 6. The same as in figure|5] but for the exchange anisotropy J^/Jx = 0.9. 



a neighborhood of the ground-state boundary between the phases |I) and |II) where a In this respect, the zero- 

temperature phase transition between the phases |I) and |II) might be regarded as a phase transition of the infinite order 
and hence, this special critical point in fact represents a quite remarkable quantum critical point. Another interesting 
finding stems from a direct comparison of the figures |5]3)-[7j3). In agreement with the aforedescribed analysis, the 
reentrant phase transitions occur in a vicinity of the zero-temperature transition between the phases |I) and |II) by 
considering the XXZ exchange anisotropy (figure [S})), while the reentrant phenomenon obviously vanishes when 
considering the sufficiently strong XYZ exchange anisotropy (figure [TJi). A disappearance of the reentrant transitions 
occurs on behalf of the XYZ exchange anisotropy, which generally lifts a macroscopical degeneracy of the phase |I), 
as it has been reasoned by the ground-state analysis. Owing to this fact, the critical lines of both four-fold degenerate 
phases |I) and |II) should always tend to zero temperature with the infinite gradient at the ground-state boundary 
between them whenever the XYZ exchange anisotropy is considered. The quite interesting situation thus appears if 
there is a small but non-zero exchange anisotropy between both interactions in the ;icy-plane (see figure |6}5). In this 




FIGURE 7. The same as in figure|5] but for the exchange anisotropy JyjJx. ~ 0.5. 



particular case, both critical lines meet at the ground-state boundary between the phases |I) and |II) with an infinite 
gradient but they curve in the same direction at small enough temperatures, which gives rise to another type of the 
reentrant phenomenon [28J. Naturally, the greater a difference between the exchange interactions Jx and Jy is, the 
smaller is a temperature range where the reentrant phenomenon may be observed. 



CONCLUSIONS 

In this work, the critical properties of the spin- 1/2 Ising-Heisenberg model with the pair XYZ Heisenberg interaction 
and two quartic Ising-type interactions have been studied in particular within the exact mapping technique based on 
the generalized star-square transformation. This algebraic transformation establishes an exact mapping equivalence 
between the proposed Ising-Heisenberg model and the Baxter's zero-field (symmetric) eight-vertex model, which has 
been subsequently used for obtaining several interesting exact results for the ground-state and finite-temperature phase 
diagrams, phase transitions and critical phenomena. The most interesting finding to emerge from the present study 
is an exact evidence of the quantum critical point, which appears whenever the weaker exchange interaction in the 
xy-plane equals the exchange interaction 7, along the quantization axis, i.e. whenever /, = inf{/r,7y}. It turns out 
that the critical exponents exhibit a peculiar singular behavior in a close vicinity of this quantum critical point, which 
consequently represents a phase transition of the infinite order 

The main emphasis of the present paper was to provide a deeper insight into how a spatial anisotropy in two quartic 
Ising-type interactions and the exchange anisotropy in the XYZ Heisenberg pair interaction influence the striking 
weak-universal critical behavior. It has been shown that a spatial anisotropy in two quartic Ising-type interactions has 
less pronounced effect on the critical behavior, which changes quantitatively rather than qualitatively upon varying the 
anisotropy in two quartic Ising-type interactions. On the other hand, the exchange anisotropy in the XYZ Heisenberg 
pair interaction significantly changes the overall critical behavior, namely, it enables to change a location of the 
quantum critical point as well as an appearance (or disappearance) of the reentrant phase transitions. 
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